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Abstract. This paper is the first of a series reporting the re- 
sults of the Australia Telescope ESO Slice Project (ATESP) 
radio survey obtained at 1400 MHz with the Australia Tele- 
scope Compact Array (ATCA) close to the South Galactic Pole 
(SGP) over the region covered by the ESO Slice Project (ESP) 
galaxy redshift survey (8 ~ —40°). The survey consists of 16 
radio mosaics with ~ 8" x 14" resolution and uniform sensi- 
tivity (la noise level ~79 /jJy) over an area of ~ 26 sq. de- 
grees. Here we present the design of the survey, we describe 
the mosaic observing technique which was used to obtain an 
optimal combination of uniform and high sensitivity over the 
whole area, and the data reduction: the problems encountered 
and the solutions adopted. 
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1. Introduction 

Sky surveys play an important role in astronomy. Large collec- 
tions of objects allow reliable studies of the average properties 
of different constituents of the universe; new populations can 
be found and well defined sub-samples can be extracted for 
further detailed analysis. 

Surveys in the radio domain have yielded many important 
results in the past, since the discovery of radio galaxies and 
quasars. The first samples of radio sources (5^1 Jy) have 
demonstrated that classical radio galaxies are rare in the local 
universe and strongly evolve with cosmic time both in density 
and luminosity (e.g. Longair 1966). More recently, deep radio 
surveys (S 1 mJy) have shown that normalized radio counts 
show a flattening below a few mJy, corresponding to a steep- 
ening in the actual observed counts (see e.g. Windhorst et al. 



199C for counts at 1.4 GHz). This change of slope is generally 



interpreted as being due to the presence of a new population of 
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radio sources (the so-called sub-mJy population) which does 



not show up at higher flux densities (see e.g. Condon 1989). To 



explain the new faint radio population several scenarios have 
been invoked: strongly-evolving normal spirals (Condon 1984 , 



1989| ) ; acti vely star forming galaxies (e.g. Rowan-Robinson 



et al. |1993|); or a non-evolving populat ion of local (z < 0.1) 



low-luminosity galaxies (e.g. Wall et al. |1986[ ). The true nature 
of the population is not well established. The same is true for 
the relative contributions from the above mentioned objects. 
Furthermore the source space density inferred from the faint 
end of the local bivariate luminosity function for both spiral s 
and ellipticals is not well known (see e.g. Condon 1996 ). 
Therefore it is not possible to estimate the local contribution 
to the counts (even if expected to be small) nor is there a 
clear local reference frame for understanding evolutionary 
phenomena. 

Unfortunately, due to the long observing times required to 
reach faint fluxes, the existing samples in the sub-mJy region 
are generally small. Table [l] shows a compilation of the largest 
1 .4 GHz surveys available in the mJy and sub-mJy regime with 
the surveyed areas and limiting fluxes (note, however, that the 
quoted limiting fluxes are often not uniform over the entire 
areas). 

The identification work and subsequent spectroscopy are very 
demanding in terms of telescope time. Typically, no more than 
~ 50 — 60% of the radio sources in sub-mJy samples have been 
identified on optical images, even though for the p]y survey 
in the Hubble Deep Field an ident ificati on rate of about 80% 
has been reached (Richards et al 



1999) 



On the other hand, 
the typical fraction of spectra available is only ~ 20%. The 
best studied sample is the Marano Field, where ^45% of the 



sources have spectral information (Gruppioni et al. 1999a ). 
To establish a firm point in the radio properties of galaxies in 
the local (z < 0.2) universe it is necessary to survey a large area 
in the sky down to faint flux limits. Furthermore it is necessary 
to have in the same region a statistically significant sample of 
galaxies with well studied optical properties (radial velocities, 
magnitudes etc.). To alleviate the identification work, regions 
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Table 1. 1.4 GHz mJy and sub-mJy radio surveys. 



Survey 


References 


Area 


Slim 






sq. deg. 


mJy 



NVSS 
ELAIS N b 
FIRST 
ELAIS S 
VLA-NEP 
PDF 

Marano Field 
ELAIS N a 
LBDS 

Lockman Hole 
Lynx 3A 
0852+17 
1300+30 
HDF 

ATESP 



Condon et al . |1998 
Ciliegi et al . |l99Sf 
White etal.flW 



Gruppioni et al. |1999b| 
Kollgaard et a l |199 
Hopkins et al. |1998| 
Gruppioni et al. |l99 
Ciliegi et al. |l999[ 



Windhorst et al. 
de Ruiter et al. 
Oort 



1987 



1984 



1997 



Condon & Mitchell 
Mitchell & Co ndon 
Richards |l999j 



1984 



1985 



this paper 



3 xlO' 
4.22 
1550 

4.0 
29.3 

3.0 
0.36 
0.12 

5.5 
0.35 

0.8 
0.32 
0.25 

0.3 

25.9 



,4 



2.5 
1.15 
1.0 
0.4 
0.3 
0.2 
0.2 
0.135 
0.1 -0.2 
0.12 
0.1 
0.08 
0.08 
0.04 

0.47 



with deep photometry (possibly multicolor) already available 
provide a significant advantage. The region we have selected 
fulfills these requirements at least partially. 
Vettolani et al. (1997) made a deep redshift survey in two 
strips of 22° x 1° and 5° x 1° near the SGP by studying 
photometrically and spectroscopically nearly all galaxies down 
to bj ~ 19.4. The survey, yielding 3342 redshifts (Vettolani 

0.1 with 10% of the 



et al. 1998 1, has a typical depth of z 
objects at z > 0.2 and is 90% complete. In the same region lies 
the ESO Imaging Survey (EIS, Nonino et al. |1999[) Patch A 



(3.2 sq. degr), consisting of deep images in the I band out of 
which a galaxy catalogue 95% complete to / = 22.5 has been 
extracted. Further V band images are available over ~ 1 .5 sq. 
degr. 

We used the 6 km configuration of the ATCA to make a 20 
cm radio continuum mosaic of the region covered by the ESP 
galaxy redshift survey. The ATESP radio survey has uniform 
sensitivity (la noise level ~79 /jJy). 

The present paper essentially deals with a description of 
the survey, the observations, the mosaic technique and the data 
reduction. It is organized as follows. In Sect. ^ the survey de- 
sign, in particular with respect to the mosaic technique is ex- 
plained. In Sects. H and Q we present in detail the calibration 
of our 20 cm observations and the data reduction. We discuss 
the problems encountered and the solutions adopted. Sect. || is 
dedicated to the analysis of the mosaics. A summary is given 
in Sect. [| 

2. Survey design 

The radio observations were carried out with two main goals 
in mind. The first aim was to detect the ESP galaxies in order 
to derive the 'local' (z ~ 0.1) bivariate luminosity function. We 



therefore tried to keep the sensitivity as uniform as possible 
over the whole ESP area, while at the same time reaching flux 
densities well below ~ 1 mJy (see Sect. 2.2). 
The second aim was to have a complete catalogue of faint ra- 
dio sources in order to study the sub-mJy population through 
a programme of optical identification of complete radio source 
samples extracted from the ATESP survey, exploiting the avail- 
able data, i.e. deep CCD images. 

As the survey is intended to achieve uniform sensitivity over 
a large area it is necessary to make use of the mosaicing tech- 
nique. 

2.1. Mosaicing Technique 

Mosaicing is the combination of regularly spaced multiple 
pointings of a radio telescope which are then linearly combined 
to produce an image larger than the radio telescope's primary 
beam. The linear mosaicing consists of a weighted average of 
the pixels in the individual pointings, with the weights deter- 
mined by the primary beam response and t he expected noise 
level (see equation in Sault & Killeen 1995 ). If the observing 
parameters and conditions are the same for every individual 
pointing, the expected noise variance in any observed field can 
be assumed to be equal for every observing field and the inten- 
sity distribution in the mosaiced final image, /(/,m), is modu- 
lated only by the primary beam response: 



I(l,m) = 



Y.jP{l-k,m- mi)lj (I , m) 
ZiP^l-lum-mi) 



(1) 



where the summation is over the set of pointing centres (li,nij). 
Ii(l,m) is the image formed from the z'-th pointing (not cor- 
rected for the primary beam response) and P(l,m) is the pri- 
mary beam pattern. 

In planning a mosaicing experiment, the main issue to be de- 
cided is the pointing grid pattern, i.e. geometry and pointing 
spacings. For a detection experiment on a large area of sky (like 
the ATESP survey) the main requirement is uniform sensitivity 
over the entire region together with high observing efficiency. 
Such requirement can be satisfied by choosing opportunely the 
pointing grid pattern. 

The mosaic noise standard deviation, a(I(l,m)), can be ob- 
tained from the error propagation of Eq. [I] and the uniform sen- 
sitivity constraint is expressed by 



a(/(/,m)) 



1 



= constant 



(2) 



Ztl?(!-lt,m-mt) 



for every position (l,m). In other words, the mosaic sensitivity, 
expressed in terms of as (a is the noise expected in the individ- 
ual pointings) is modulated by the squared sum of the primary 
beam response. 

The ATCA primary beam pattern can be approximated by a 
circular Gaussian function (Wieringa & Kesteven 1992) 



P,W=e" 41n HFWHPy 



(3) 
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where r = yP- + m 2 is the radial distance from the image phase 
center and FWHP is the full width at half power of the pri- 
mary beam (~ 33 arcmin for ATCA observations at 1.4 GHz). 
Since the square of a Gaussian is still a Gaussian with FWHP 
reduced by a factor \f2, the quadratic sum in Eq. ^ can be writ- 
ten as a linear sum of Gaussians with FWHP' = FWHP/\/2 
(corresponding to s ~ 33'/V5 = 23.3' for ATCA 1.4 GHz ob- 
servations). 

To make the final choice for the ATESP survey grid pattern, we 
have performed a series of simulations of the bidimensional 
quantity 

V^KfM*) 

where s t = r = (I - U) 2 + (m - im) 2 and FWHP' = 

23.3', varying the pointing spacings s and the grid geometry 
(hexagonal and/or rectangular grids). 

In general a very good compromise between uniform sensitiv- 
ity and observing efficiency is represented by a grid pattern 
with pointing spacings of the order of s ~ FWHP'. For our par- 
ticular case, the best choice turned out to be a 20' spacing rect- 
angular grid. The mosaic noise variations over a reference area 
of 1 sq. degr. for such a grid configuration are shown in Fig. Ill 
As expected, in the region of interest (central box) the noise 
is rather constant: variations are < 5%, except at the region 
borders (< 10%). We point out that hexagonal grids should be 
preferred when imaging wide areas of sky (like for the NVSS 
and FIRST radio surveys), but are not very efficient in the case 
of a narrow 1 -degree wide strip of sky. 

2.2. Observing Frequency, Resolution and Sensitivity 

At 20 cm (the longest observing wavelength available at the 
ATCA) the field of view is largest (gaussian primary beam 
FWHP~ 33') and the system noise is lowest. Thus, observ- 
ing at 20 cm allows minimization of both the number of fields 
(i.e. pointings) required to complete the survey and the observ- 
ing time spent on each field. The ATCA can observe at two 
frequencies simultaneously (for instance 20 and 13 cm). How- 
ever, we decided to optimize sensitivity at the expense of spec- 
tral information, by setting both receivers in the 20 cm band 
and observing in continuum mode (2 x 128 MHz bandwidth, 
each divided into 32 x 4 MHz channels in order to reduce the 
bandwidth smearing effect). This choice was also influenced by 
another consideration: since the field of view depends on the 
observing frequency, the grid pattern could not be optimized to 
get uniform sensitivity for both 13 and 20 cm bands simultane- 
ously. 

The observations were carried out at full resolution (ATCA 6 
km configuration), since the identification follow up benefits 
from high spatial resolution, and the expected fraction of very 
extended sources (that could be resolved out and lost) is low at 
the ATESP resolution. Using the angular size distribution given 
by Windhorst et al. (1990) for radio sources, we estimate that 




RA offset (pixels) 



Fig. 1. Noise variations expected over a 1 sq.degr. mosaic (cen- 
tral box, 1 pixel = 1') for a 20' spacing rectangular grid. Filled 
circles indicate the pointing centres. At the center of the mo- 
saic o(7) /o ~ 1 . Contours refer to 5, 10, 20, 40, 80, 160, 320% 
noise increment from the center. 

70% of the mJy and sub-mJy sources would appear point- 
like at the ATESP resolution, and 5% would have angular 
sizes twice the beam size or larger. 

Since the ESP sample distance distribution peaks at z — 0.1 
and 'normal' galaxies are typically low-power radio sources, 
deep radio observations were needed to ensure detections of 
a statistically significant number of ESP galaxies. We con- 
sidered satisfactory a point source radio limit of the order of 
~ 0.2 mJy (3a), which corresponds to a detection threshold of 
P ~ 3 x 10 21 W Hz" 1 at z ~ 0.1 (H G = 100 km/s Mpc" 1 ). Fur- 
thermore a large sample of sub-mJy radio sources can be con- 
structed at a 6o detection limit, corresponding to a flux limit of 
- 0.5 mJy. 

3. Observations and calibration 

To cover the two areas of 22° x 1° (region A) and 5° x 1° (re- 
gion B) of the ESP survey with uniform sensitivity, 69 x 4 and 
15x4 pointings at 20' spacing are needed respectively. An area 
of 1.3 sq. degs. (4x4 pointings) of region A was not observed, 
because of the presence of the strong radio source PMNJ2326- 
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Table 2. Log of the observations. 



4.1. Flagging of bad data 



Date 


f obs 


Array 


Vi 


v 2 




h 




MHz 


MHz 


18/11/94-21/11/94 


3 x 12 


6D 


1344 


1452 


23/12/94-04/01/95 


13 x 12 


6A 


1344 


1452 


15/12/95-01/01/96 


17 x 12 


6C 


1344 


1448 



4027, which would have prevented us from reaching the deep 
noise level required. This reduced the total number of fields to 
be observed to 65 x 4 (region A) and 15x4 (region B), i.e. 320 
pointings in total. The right ascension range actually covered 
by the ATESP survey is thus 23 h 31 m - 01 h 23 m and 22 h 32 m to 
22 h 57 m (J2000). 

An observing time of 1.2 hours per pointing was needed to 
reach a 3o RMS noise level of ~ 0.2 mJy (with 2 x 128 MHz 
bandwidth). To obtain good hour angle coverage we organized 
the snapshot observations as follows: the 320 fields were di- 
vided in 16 sets of 20 fields each (5 x 4). Every set was then 
observed for 2 x 12 h , that is, the 20 fields were observed in 
sequence, for 1 minute each, repeating this for 72 times and 
adding 3 minutes for calibration every hour. 
The observing campaign started in November 1994 and was 
completed in January 1996. A log of the observations is given 
in Table || (dates, arrays used, observing time and frequencies). 
We stress that, for our purposes, the use of different 6 km ar- 
rays is not relevant in any way. 

The two 128 MHz observing bands were set in the most 
interference-free region of the 20 cm band (1.3—1.5 GHz). 
The flux density calibration was performed through observa- 
tions of the source PKS B1934-638, which is the standard pri- 
mary calibrator for ATCA observations (S = 14.9 Jy at v = 
1384 MHz, Baars et al. 1977 flux scale). The phase and gain 
calibration was based on observations of secondary calibrators, 



selected from the ATCA calibrator list. Every single 12 run 
and each of the two observing bands were calibrated separately 
following the standard procedures for ATCA observations. 



4. Data reduction 

For the data reduction we used the Australia Telescope Na- 
tional Facility (ATNF) release of the Multichannel Image Re- 
construction, Image Analy sis and Display (MIRIAD) software 
package (Sault et al. 1995 ). A number of steps in the reduction 
process are usually done interactively, specifically the removal 
of bad data ('flagging'), and, in a certain measure, cleaning and 
self calibration. Due to the large amount of data involved we 
found it more practical to develop a semi-automated reduction 
pipeline. 



We made a modified version of the MIRIAD task TVFLAG (in- 
serted in MIRIAD as TVCLIP), which recursively flags visibil- 
ities with amplitudes exceeding a given threshold. The thresh- 
old was set as a convenient multiple of the average absolute 
deviation (|A5|) from the running median, evaluated separately 
for each baseline, each channel and each integration cycle (10 
s). 

For the primary calibrator the automated flagging procedure 
was applied before the calibration. This was necessary to avoid 
the calibration being affected by bad data. For the secondary 
calibrators and the mosaic data the bandpass and instrumental 
polarization calibration were applied before running TVCLIP. 
As we noticed that the shortest baselines introduced some low 
level, spatially correlated features in the images, which could 
affect the zero level for faint sources, we decided to remove all 
baselines shorter than 500 m from the data prior to the pipeline 
processing (rejection of ~ 10% of the visibilities). As a conse- 
quence, the ATESP survey becomes progressively insensitive 
to sources larger than 30": assuming a Gaussian shape, only 
50% of the flux for a 30" large source would appear in the 
ATESP images. However, the expected fraction of sources with 
angular sizes > 30" is very small: < 2% at fluxes S < 1 mJy ac- 
cording to the Windhorst et al. (1990) angular size distribution. 

4.2. Cleaning and Self-calibration 

Since the primary beam response is frequency dependent, we 
did not merge the data from the two observing bands before 
imaging and cleaning. This results in a slightly poorer UV cov- 
erage but allows the cleaning process to succeed in subtracting 
correctly 100% of the source flux, and self-calibration to be 
more effective and reliable. 

On the other hand, to improve UV coverage and sensitivity, for 
each field we have merged the (calibrated) data coming from 
all the observing runs. 

In contrast with the imaging of extended sources, joint de- 
convolution is not needed for a point source survey. It is also 
very expensive computationally for high resolution images. We 
therefore reduced every field separately, simplifying imaging 
considerably. 

For each field a 2048 x 2048 pixel image (total emission only) 
was produced (pixel size = 2.5"). The entire image was then 
cleaned in order to deconvolve all the sources in the field. To 
improve sensitivity we used natural weighting, which gave a 
synthesized beam typically of the order of 8" x 14". 
Each image went through different cleaning cycles. First, we 
produced the list of the brightest components to use as model 
for self-calibration. Phase only self-calibration was applied. 
Usually two iterations were sufficient to remove phase error 
'stripes' and improve the image quality. 

The self-calibrated visibilities were then used to produce a 
deeper cleaned image. A serious problem arises if snapshot im- 
ages are cleaned too deeply. Due to the incomplete UV cov- 
erage the number of CLEAN components can approach the 
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Table 3. Main parameters for the 16 final mosaics. 



Mosaic" 


Fields 


Tangent Point* 


Synthesized B 


eam e 






< rj > 


fid x to y 


nx m 




on 


^min * ^maj 


PA. 


mJy 


<"Jy 


<"Jy 


fld01to06 


6x4 


22 35 57 37 


-39 SO 1 S 


7.8" X 


12.9" 


1° 


-0.49 


78.7 


78.5±3.0 


fld05toll 


7x4 


22 44 57.60 


-39 59 15 


7.8" X 


12.8" 


3° 


-0.96 


77.8 


78.5±3.8 


fidl0tol5 


6x4 


22 53 21.68 


-39 59 15.0 


7.9" x 


13.0" 


-1° 


-0.70 


88.1 


84.1±5.3 


fld20to25 


6x4 


23 34 56.09 


-39 58 48.0 


8 5" * 

O.J A 


16 8" 

1U.O 


1° 


-0.68 


83.0 


80.3±4.3 


fld24to30 


7x4 


?^ 38 76 


-39 58 48 

J7 JO T^O.U 


9.5" x 


16.6" 


-2° 


-0.50 


82.8 


83.2±6.4 


fld29to35 


7x4 


23 52 20.43 


-39 58 48.0 


8.9" x 


16.7" 


2° 


-0.53 


79.2 


76.2±2.9 


fld34to40 


7x4 


00 01 02.59 


-39 58 48.0 


8.0" x 


14.5" 


5° 


-0.48 


76.3 


74.8±4.0 


fld39to45 


7x4 


00 09 44.76 


-39 58 48.0 


8.6" x 


14.3" 


-10° 


-0.61 


81.2 


80.4±3.4 


fld44to50 


7x4 


00 18 26.93 


-39 58 48.0 


7.4" x 


14.0" 


7° 


-0.43 


78.0 


76.4±1.7 


fld49to55 


7x4 


00 27 09.10 


-39 58 48.0 


8.0" x 


12.9" 


10° 


-0.62 


78.6 


77.1 ± 1.4 


fld54to60 


7x4 


00 35 51.26 


-39 58 48.0 


8.2" x 


12.6" 


4° 


-0.42 


77.3 


75.8± 1.3 


fld59to65 


7x4 


00 44 33.42 


-39 58 48.0 


7.7" x 


12.5" 


4° 


-0.44 


79.4 


77. 8± 1.4 


fld64to70 


7x4 


00 53 15.58 


-39 58 48.0 


7.7" x 


12.6" 


1° 


-0.45 


75.1 


74.6±2.2 


fld69to75 d 


7x4 


01 01 57.70 


-39 58 48.0 


7.5" x 


12.8" 


-2° 


-0.47 


81.9 


79.3±2.1 


fld74to80 


7x4 


01 10 39.86 


-39 58 48.0 


7.0" x 


15.0" 


1° 


-0.61 


77.1 


76.3±3.8 


fld79to84 


6x4 


01 19 22.03 


-39 58 48.0 


6.8" x 


14.3" 


7° 


-0.41 


68.9 


67.6±3.5 



" x and y refer to the first and last field columns composing the mosaic. J2000 reference frame. 

c PA. is defined from North through East. " Reported values for S m j n and noise refer to masked mosaic (see text). 



number of independent UV points. At this point the cleaning 
algorithm is not well constrained anymore and can interprete 
noise (sidelobes, calibration errors, etc.) as CLEAN compo- 
nents. This process can redistribute the noise into the sources 
and, as a result, the process does not converge and, in princi- 
ple, the image can be cleaned to zero flux. This produces many 
faint spurious sources, while the flux of real sources is sys- 
tematically un derest imated. This effec t has b een mentioned by 
Condon et al. ( |1998| ) and White et al. ( |1997| ) for deep snapshot 
observations with the VLA and is referred to as 'clean bias'. 
From our tests we found that cleaning down to 3a, the noise 
gets ~ 20% lower than theoretical. Down to 2a it is a factor of 
two lower and at la can be 4-5 times lower. Thus we decided 
to stop any cleaning process when the peak flux residuals are 
of the order of 4-5 times the theoretical value (setting a cut-off 
of 0.5 mJy) to minimize the clean bias (a few percent effect ex- 
pected, but see discussion in Sect. 5.3). 

After this first phase of self-calibration and deep cleaning, we 
subtracted the sources from the visibility data and we repeated 
the flagging procedure on the residual visibilities in order to re- 
duce the 'birdies' effect mentioned in Sect. 



5.1 



below. 



We then proceeded with another phase of cleaning, and pro- 
duced a half-resolution residual image covering 4 times the 
original area; only external parts of this image were cleaned 
(not the inner quarter, corresponding to the original field, which 
was already cleaned down to 0.5 mJy). This procedure allowed 
us to remove the sidelobes from more distant sources (belong- 
ing to adjacent fields). These new components were subtracted 
from visibility data before restoring the sources in the final 
cleaned image. 

As a final step we checked for bright, extended sources in the 



field, which needed deeper cleaning. A small box containing 
such a source was cleaned, to a 2a level. In general, the appli- 
cation of all these cleaning steps produced good quality single 
field images. 

4.3. Mosaicing 

The cleaned single field images were co-added in mosaics in 
order to improve the signal to noise ratio and get uniform sen- 
sitivity. Each set of 5 x 4 fields observed in 2 x 12 h observing 
blocks produces a separate independent mosaic; an overlap be- 
tween adjacent mosaics was created by adding one (or two) 
extra column(s) of fields to the side(s) of each mosaic. 
Before any mosaic is produced, every field was restored using 
the same values for the beam parameters. The restoring param- 
eters were chosen as the average value of all the fields compos- 
ing the mosaic at both frequencies. 

The final mosaics were obtained in two steps. First a single 
frequency mosaiced image was produced for each of the two 
observing bands, then the final mosaic was obtained by averag- 
ing (pixel by pixel) the two initial mosaics in order to improve 
the sensitivity. 

One of our final mosaics is shown in Fig. ^ as an exam- 
ple. The rectangular box indicates the region corresponding 
to the ESP redshift survey, where the radio survey was de- 
signed to give uniform noise. Such a region covers an area 
of 1.67° x 1° or 2° x 1° for mosaics composed by 6 x 4 
or 7 x 4 fields respectively (see Table ^|). Hereafter we will 
always refer to the central box only in our mosaic analy- 
sis. All mosaics are available thro ugh the ATESP page at 
http : / / www .ira.bo.cnr.it. 
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Fig. 2. One of the 6 x 4 fields mosaiced images. The rectangular box indicates the region corresponding to the ESP redshift 
survey. 
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Table 4. clean bias average corrections (see text). 




xlO- 4 



Fig. 3. Histogram of the (residual) flux in one of the 6x4 fields 
mosaiced images (fld01to06). As expected the flux is peaked at 
zero and the distribution is gaussian. 

5. Mosaic analysis 

Table |3| summarizes the main parameters for the final 16 mo- 
saics: for each mosaic are listed the number of fields compos- 
ing it (columns x rows), the tangent point (sky position used 
for geometry calculations) and the synthesized beam (size and 
position angle). The spatial resolution can vary from mosaic to 
mosaic depending on the particular array (6A, 6C or 6D) used 
in the observations. The mean value for the synthesized beam 
is ~ 8" x 14". 



5.7. Noise 

The last three columns of Table || show the results of the noise 
analysis. For each mosaic we report the minimum (negative) 
flux (S m i n ) recorded on the image (typically |5 m ; n | is of the or- 
der of 0.5 mJy, corresponding to the value at which we have 
stopped the cleaning) and the noise level. This has been eval- 
uated either as the FWHM of the gaussian fit to the flux dis- 
tribution of the pixels (in the range ±5 m i n ), in order to check 
for correlated noise ((Jfit), or as the standard deviation of the 
average flux in several source-free sub-regions of the mosaics, 
in order to verify uniformity (< o >). 

As expected, the noise distribution is fairly uniform within each 
mosaic and from mosaic to mosaic (< 10% variations). Also, 
for each mosaic, the two noise values are consistent, that is the 
noise can be considered gaussian (see also Fig. ||). 
On average the noise level is ~ 79 /jJy. The typical detec- 
tion limits for the ATESP survey are thus ~ 0.24 mJy at 3o 
and ~ 0.47 mJy at 6a. Dynamic range problems can cause 
slightly higher noise levels of ~ 100 /jJy around strong sources 
(Speak > 50 — 100 mJy). Such problems appear to be serious 



Mosaic 



fld34to40 1616 0.09 ±0.04 0.85 ±0.08 
fld44to50 2377 0.13 ±0.03 0.75 ±0.06 
fld69to75 3119 0.16±0.05 0.67±0.09 



only in one mosaic (fld69to75): the region around the bright 
radio source PMNJ0 104-3950 suffers from a very high noise 
level and a number of spurious sources are present. This was 
due to strong phase instabilities during the observations which 
could not be removed by self-calibration. This region (of size 
~ 20' x 25') was masked and therefore excluded from further 
analysis. Excluding this region, the total unifom sensitivity area 
covered by the ATESP survey is 25.9 sq. degr. or 7.9 10~ 3 sr. 



5.2. Artefacts 

Another problem we faced was the possible presence of arte- 
facts in the images, like spurious sources ('ghosts') at a level of 
0.1-0.5% opposite to bright sources with respect to the phase 
center of the image and 'holes' in the centre of the field. The 
first problem, caused by the Gibbs phenomenon, arises from 
the use of an XF correlator and can be serious in high dynamic 
range continuum observations (like ours). The second effect 
is a system error produced by the harmonics of the 128 MHz 
sampler clock at 1408 MHz. Both effects can be completely 
removed as long as the observing bands are centered appropri- 
ately (Killeen |1995[ Sault |l995| ). Unfortunately, at the time of 
our first two observing runs these effects were not yet known. 
We therefore could apply the corrections only to the data taken 
in the last observing run. 

We point out that the corrections, when applied, result in a 
larger bandwidth smearing effect, since only 13 x 8(10) MHz 
channels are used (instead of 32 x 4 MHz). 
Wherever not corrected for, the 'ghosts' problem is unlikely to 
be serious in our case, since 'ghosts' appear in different places 
for each field and so they tend to average out when mosaicing 
the fields. Moreover, only radio sources brighter than ^100 
mJy can produce detectable 'ghosts' in the final mosaics and 
such bright sources are very few in the region surveyed (~ 30) 
and therefore could be easily checked. No evident 'ghosts' have 
been found. 

We have also tried to reduce the sampler clock self-interference 
effect as far as possible by flagging residual bad visibilities 
(that is correlated noise) after a first step of cleaning and self- 
calibration (see also Sect. \.2 ), but some of our fields still show 
it to a small extent. On the other hand the area of sky affected 
by 'holes' is of the order of a few percents (< 3%) of the total 
region observed. 
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5.3. Clean Bias 




Fig. 4. The source flux measured after the cleaning (Soutput) nor ~ 
malized to the true source flux (Sj nput ) as a function of the flux 
itself (expressed in terms of c) is shown for three different 
cases: 1616 cc's mosaic (top panel), 2377 cc's mosaic (mid- 
dle panel) and 3119 cc's mosaic (bottom panel). Also shown 
are the linear fits (see text). 



2500 

cc's 



Fig. 5. We present Soutput/Sinput as a function of the average 
number of clean components. Each dotted line refers to a dif- 
ferent source flux (100a, 50a, 30a, etc.). Also shown is the av- 
erage number of clean components for each of our 16 mosaics 
(full lines at the bottom of the figure). 



As already mentioned in Sect. 4.2, when the UV coverage 
is incomplete, the cleaning process can 'steal' flux from real 
sources and redistribute it on top of noise peaks producing spu- 
rious ones. 

To quantify the actual effect in our mosaics we performed a 
set of simulations by injecting point sources in the survey UV 
data at random positions. Then the whole cleaning process was 
started. The number of sources injected in a mosaic (~ 500) 
was chosen in order to reduce the statistical uncertainty without 
changing significantly the components/image statistics. The 
source fluxes cover the entire range of the survey: 250 sources 
at 3a, 100 at 5a, 50 at 7a, 50 at 10a, 25 at 20a, 10 at 30a, 10 
at 50a, 10 at 100a. 

Taking into account the time and frequency resolution and 
the baseline lenghts, we get about 1500-2500 indepentent UV 
points for each field observed. This means that with about 2000 
independent (not too close together or on top of each other) 
clean components we could clean the image to zero flux. The 
average number of (not independent) clean components per 
mosaic ranges between 1500 and 3200. Since we used a clean 
loop gain factor g = 0. 1 , the number of independent clean com- 
ponents per mosaic can be estimated as 1/10 of the numbers 
reported above. We thus expect a significant clean bias effect 
(10-20%), larger for mosaics with a higher number of cleaning 
components. We then decided to test three mosaics, with a low, 
an intermediate and a high average number of cleaning com- 
ponents (cc's) respectively: fld34to40 (1616 cc's), fld44to50 
(2377 cc's) and fld69to75 (3119 cc's). 

The results of the tests are presented in Fig. Q and Fig. || and 
summarized in Table Q Fig. ^ shows, for each of the three 
mosaics, the average source flux measured after the cleaning 
(^output) normalized to the true source flux (5i nput ) as a func- 
tion of the flux itself (expressed in terms of a). In general the 
clean bias increases going to fainter fluxes and, as expected, 
depends on the number of cleaning components. In the best 
case (1616 cc's) we get < 10% flux underestimation for the 
faintest sources; in the worst case (3119 cc's) the effect rises 
up to ~ 20%. 

The dependency of the clean bias on the number of clean- 
ing components is more clearly shown in Fig. |[ Here we 
present Soutput/Sinput as a function of the average number of 
clean components for different source fluxes (100a, 50a, 30a, 
etc.). Again, the clean bias affects more seriously the faintest 
sources. Moreover it is evident that we are not dealing with a 
linear effect: a sudden worsening appears at fluxes of the order 
of 10-20C and when the number of cleaning components ex- 
ceeds - 2000. 

A first order fit of the clean bias effect for the three differ- 
ent mosaics has been obtained by applying the least squares 
method to the function Soutput/Sinput = fl l°g(^input/o) +b. The 
values obtained for the parameters a and b are listed in Table Q 
and the curves are shown in Fig. 0. 
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6. Summary 

The ATESP survey at 1.4 GHz is based on snapshot obser- 
vations of 320 overlapping primary beam fields, reduced 
separately and then combined together to produce 16 big mo- 
saiced images. The total area surveyed with uniform sensitivity 
(la noise level ~ 79/jJy) covers 25.9 sq. degr. The spatial 
resolution is typically <~ 8" x 14" providing radio positions 
with internal accuracy of the order of 1 arcsec for 6o radio 
sources. 

We stress the importance of estimating the relevance and the 
behaviour of the so called clean bias effect for any deep survey 
obtained with snapshot observations. The source fluxes can be 
seriously affected by such problem and for reliable scientific 
analysis the effect must be taken into account and corrected for. 

Future papers in this series will deal with: 

1 . the radio sources catalogue complete down to a limiting 6o 
flux density of ~ 0.5 mJy 

2. the ATESP radio source counts 

3. the radio properties of the ESP galaxies and the local bi- 
variate luminosity function 

4. the optical identifications and spectroscopy of the objects 
in the EIS sub-region. 
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